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Abstract 

Gyrosonic is a novel audio binaural stimulus, that generates rotational perceptions of 
sound movement in the brain at a particular predetermined frequency The influence of 
gyrosonic on the autonomic nervous system of healthy subjects has been examined by 
analyzing heart rate variability (HRV) in time- and frequency- domain. The M-lagged 
Poincare plot shows that the parameters SDl, SD2 and ratio SD12 (SD1/SD2) increases 
with lagged number M, and M-dependence is well described by Pade' approximant X ij^^^f 
where values of X) /3 and 7 depends on parameters SDl, SD2 and SD12. The values of 
parameters SDl, SD2 and SD12 after gyrosonic stimulation are augmented, compared to 
pre-stimulation values for all M. The slope and curvature that deflne the variation of SD12 
with M increase considerably due to the stimulation. Slowing down of the Heart rate and 
increase in the standard deviation SD and the root-mean squared successive differences 
(RMSSD) in the RR-interval are observed after stimulation. The spectral power of both 
low (LF) - and high (HF) - frequency go up due to the stimulation. The DFA analysis of 
RR interval exhibits a decrease in value of exponent a due to the gyrosonic. The results 
indicate that there may be an improvement of the sympatho-vagal balance due to this 
novel audio stimulus. 
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1 Introduction 



Autonomic nervous system (ANS), with its main two divisions; Sympathetic and Parasym- 
pathetic, may be viewed as a hierarchically coordinated neuronal network. It regulates the 
exchange of energy and information with the environment [1]. Most of the internal organs like 
the heart, gastrointestinal tract, lungs, urinary bladder and blood vessels are under the influ- 
ence of the autonomic nervous system. Usually the sympathetic and parasympathetic nervous 
systems have opposing effects on the organs. The degree of influence may vary and this deter- 
mines the ultimate activity of that organ. The sympathetic nervous system is rapidly activated 
in physical or mental stressful conditions. It increases the heart rate, cardiac output and blood 
flow to the muscles. It dilates the pupil and decreases gastrointestinal tract activity. Its action 
is therefore, sometimes aptly referred to as the 'flight or fight' response. On the other hand, 
the parasympathetic system decreases the heart rate and blood pressure, constricts the pupil 
and increases gastrointestinal activity. The 'rest and digest' response is coined for its action [2]. 
The body maintains a proper balance between the sympathetic and parasympathetic activity 
for appropriate functioning. A deviation from this balance may result in disease conditions like 
acute coronary syndrome, chronic heart failure, diabetes mellitus etc [3]. 

It has been shown that sensory inputs, like a harmonic auditory stimulus, can have a wide 
range of psychological and physical effects [4]. A sensory input can trigger a cascading series of 
events along the nerve down which it travels [5]. Many of the beneficial effects of the auditory 
stimulus take their origin along the route of the impulse. Music has the capacity to mod- 
ify the psychobiological state of human. It can thus relieve stress and stress related ailments 
[6].Gyrosonic is a novel auditory stimulus and appears to possess such capacity. When applied 
binaurally, it produces the perception of rotation of the audio source inside the brain. The 
frequency of this rotation is in the infrasonic region (~ 2Hz). The perception of movement in 
the auditory space by human depends on a number of cues. It has been demonstrated that the 
moving sound produced by the sequential excitation of a mono source through several speakers 
in a free field, can generate a specific activity in the brain [7] . Sound motion has been shown to 
evoke a magnetic field inside the brain, and the evoked magnetic response is specific to moving 
sound and is absent for non- moving sound stimulus [8] . Another study has shown that the right 
parietal cortex is involved in the processing of sound motion [9]. Compared to the earlier ex- 
periments gyrosonic stimulus has some improved rotational features . Thus it can be expected 
that this may produce a larger spectrum of brain activation. Earlier studies with this moving 
sound showed that the arousal level of psychosomatic patients was significantly reduced [10]. 

The heart rate variabihty (HRV), which is a measure of the beat to beat fluctuation of the 
heart rate, reflects the time varying influence of ANS and its components on the cardiovascular 



function. Due to its non-invasive nature and convenience of the measurement, the HRV is often 
used to assess the influence of the autonomic nervous system on the cardiac function. Under 
normal circumstances the HRV is regulated by the combined action of heart automaticity and 
the feedback elements from the vagal and sympathetic activity on the heart. The stimuli that 
are capable of altering the feedback components therefore change the dynamics of the heart. 
Gyrosonic stimulus is assumed to be one such stimulus. In this paper we have estimated the 
influence of gyrosonic stimulus on the dynamics of the heart by assessing the Heart Rate vari- 
ability (HRV) in response to this stimulus. The HRV can be assessed in time domain and in 
frequency domain. The time domain analysis provides quantitative measurement of the varia- 
tion of heart rate, standard deviation SD, root mean square successive differences (RMSSD) in 
RR-interval.On the other hand the frequency domain analysis of HRV mainly measures the high 
frequency (HF) and low frequency (LF) power spectrum of fluctuation of the RR interval. It 
has been suggested that the high frequency power spectrum is modulated by parasympathetic 
(mainly vagal) activity and the low frequency power spectrum is mainly influenced by both 
parts of ANS activity [11,12]. 

Dynamics of the cardiac system is nonlinear. Hence it is reasonable to assume that a 
nonlinear analysis is more appropriate means to get an accurate idea about the cardiac system. 
Poincare' plot analysis meets this criterion [13, 14, 15]. It is a nonlinear geometric method of 
heart rate analysis. It is basically a scatter plot of any heart rate interval RRn and next one 
i?-RAr+i.When this plot is adjusted with an ellipse, three important parameters then define the 
plot. These are the length of the semi-minor axis of the ellipse SDl, the length of the semi-major 
axis of the ellipse SD2, and their ratio SD12.The parameter SDl is the standard deviation of 
instantaneous beat to beat heart rate variability and it is the measure of short term variability 
of heart. It is mainly influenced by parasympathetic regulation on the heart. Other parameter 
SD2 is the measure of long term variability. It has been shown that these parameters are 
correlated with the power spectral density of heart rate fluctuation. SDl is correlated with 
the high frequency spectrum and SD2 is mainly associated with the low frequency spectrum 
[16]. Instead of taking the successive RR interval, the Poincare'plot can be generalized by taking 
a lag of greater than 1 between RR intervals. A study with a lag of 4 has shown evidence of 
being superior to normal Poincare' [17]. Moreover, a study on chronic heart failure patients has 
demonstrated strong dependency of the parameter SD12 on the lag interval. This study has also 
shown that the curvature of the plot of SDl 2 with various lag numbers is signiflcantly different 
in patients and normal subjects. In chronic heart failure patients the curvature is much smaller 
[18]. Another study on diabetic patients has pointed out that the SDl decreases in diabetic 
patients with higher lag numbers [19]. Using this lagged Poincare' plot it has been reported 
that SDl and SD2 decreases with an increase in the lag number for smokers [20]. Apart from 



the spectral analysis and Poincare' plot analysis for the heart rate variability, another method, 
called Detrended Fluctuation Analysis, can be applied with some added advantages. It has been 
found that under healthy conditions, the RR interval time series exhibits long-range power law 
correlation. Similar behavior has also been found in other physical systems [21]. The Detrended 
Fluctuation analysis proposed by Peng et al is supposed to reflect on this long range correlation 
of the RR interval time series [22, 2 3] and is widely used. 

In our study we have used all these methods to examine the change in HRV in healthy 
subjects after the gyrosonic stimulation. Through this study we propose to show that a change 
in the autonomic regulation occurs due to the influence of the gyrosonic stimuli. 

2 Stimuli 

In this study we have used a gyrosonic stimuli constructed from an Indian percussion instru- 
ment, the Tabla. The stimuli were recorded digitally at a sampling rate of 44.1 kHz and in 16 
bit. Then the amplitude of sound was modulated at 2 Hz in such a way that the amplitude 
increased in one ear and decreased in the other ear. This modulation created a perception of 
moving sound. When one ear heard the sound with higher amplitude it assumed the sound 
to be moving towards that ear whereas the ear hearing the lower amplitude sound took that 
the sound was moving away. This produced the effect as though the sound was moving in a 
horizontal plane around the head. There was a phase of 0.568 seconds between the advancing 
or the rising and receding or the falling stimuli. The slopes of the amplitude of the rising or 
falling sound were exponential. The stimulus was presented through headphones, connected to 
a computer through a custom electronic interface. The sound was played at a sensation level 
around 50 db. The ultimate perception produced was that of a rotating sound moving in a 
horizontal plane inside the brain. The duration of the stimuli was 9.5 minutes. This duration 
was chosen to ensure that the subject perceived the effect in its optimum potential. No subject 
experienced nausea or episode of vertigo during or after the playback of the stimuli. The stimuli 
were given in a sound proof room with the subjects in recumbent position. 

3 Subjects and Measurements 

Thirty one subjects (11 male, 20 female with average age 36 ± 12 yrs) volunteered for the 
study. All of the subjects were explained about the study and all of them gave their consent 
for the study. All the subjects were healthy and were not on any kind of medication. There 
was no history of any type of nervous system disorder. The study was approved by the ethical 
committee of IIT Kharagpur, India. ECG data was taken from three limb leads. Patients 



were in supine condition. Sampling rate was 500 Hz with a resolution of 12 bit. A pre stimuli 
ECG was taken for 10 mins. Then the subjects hstened to the gyrosonic stimulus. The post 
stimulation ECG for 10 mins was recorded after a gap of 10 mins from the end of the stimulation. 
Whole 10 mins data of ECG for both pre and post stimuli was analyzed after selecting the sinus 
beats only. Ectopic beats were detected visually and deleted manually. RR interval or better 
to say normal beat to normal beat interval was detected through the Origin software. The 
non linear analysis like the Poincare' plot and the Dretended Fluctuation analysis were done 
with program written in Matlab Software. The linear HRV analysis was done using standard 
program. 

4 Data analysis 

During the analysis of HRV wc have strictly followed the Task Force guideline of HRV anal- 
ysis[13]. In linear time domain the heart rate (HR), standard deviation (SD) of RR interval 
(basically the normal beat to normal beat interval), root-mean squared successive differences 
(RMSDD) were calculated. The ratio of SD/RMSDD was considered to be a good measure of 
the sympatho- vagal balance [14]. In the frequency domain, the power spectral density for the 
low frequency (0.04-0.15 Hz) and the power spectral density for the high frequency (0.015-0.5 
Hz) were calculated. The distribution of power and the central frequency of LF and HF were 
not fixed. They usually vary in relation to changes in the autonomic modulation of the heart. 
According to the Task Force Guidehne the RMSDD represented the HF variabihty and the to- 
tal variability of the RR interval data was reflected by the SD. The parameters SDl and SD2, 
as described earher, were the width and the length of the ellipse. The ellipse was fltted to the 
Poincare plot of RRN-I-M vs RRN where M is the lag number. SDl and SD2 were calculated 
for lag M from the relation :SD1 = ($(M) - $(0))^/^ and SD2 = ($(M) + $(0))V2 where the 
auto-covariance function $(M) is given by 

$(M) =E[{RRn - RR){RRn+m - RR)]- 

[16]. 

The long term correlation in RR-time sequence was assessed by Detrended Fluctuation 
Analysis [22]. The measure of correlation was given by a scahng exponent (a) of the fluctuation 
function F{t) ^ r". The computation of fluctuation function F{t) was done in the following 
way. For a given time sequence R{ti), U — iSt where St is characteristic time interval for the 
sequence and i=l,N , an integrated time series r{ti) was deflned as r{ti) — J2][R{tj) — Rm]: 
i — 1,N where R^ was the mean of R{ti). The integrated series was divided into boxes of equal 
size of time r = n5t and linear function was used to flt box data. The fluctuation function 



F(t) was calculated as root mean square fluctuations relative to the linear trend. The power 
law behaviour of F{t) provided the scaling exponent. It has been observed that acceptable 
estimate of the scaling exponent a (from DFA) can be obtained from analysis of data sets of 
length 256 samples or greater (equivalent to approximately 3.5 min for RR data at a heart rate 
of 70 bpm). The analysis of RR data for period of 10 min time interval was therefore expected 
to provide an adequate measure of the scaling exponent. 

5 Results and discussions 

The results of linear analysis of HRV is summarized in Figs.l. The Fig.l depicts the mean value 
of HR , the mean of SD and RMSSD in RR-interval of all subjects. In all figures presented here 
the term PRE and PST indicate respectively before and after gyrosonic stimulation. The heart 
rate was changed from 83.1 to 81.8 per min. The slowing down of HR was found to be more for 
subject with higher initial HR. Both SD and RMSSD were increased in post stimulation. The 
LP and HP components of spectral power in RR-interval were also increased after stimulation 
(Pig.l). The increase in LP- power was more than that of HP component. However, the ratio 
of LP/HP power was decreased in post condition. The changes of all the above parameters are 
with p < 0.05 

In Pig.2 the lagged Poincare'plot of a subject was presented with lag of 1,5,9 and 13 plot. The 
left and right part of the figure represent respectively the situation before and after gyrosonic 
stimulus. As lag number increased the plot became more scattered with consequent increase 
in both width and length of the plot. After stimulation the RRn+m vs RRn plot were more 
scattered, and center of the plot was shifted to higher value indicating slower heart rate. 

The group average values of parameters SDl, SD2 (both in sec.) and SD12 (SD1/SD2) ob- 
tained from corresponding values of individual subject were plotted against lagged number M 
(Pig. 3 points). Both parameters SDl and SD2 were increasing function of lagged number. After 
gyrosonic stimulation (PST) their values were higher than those before stimulation (PRE) and 
the growth rate of SDl with M is also higher. The result points out that the gyrosonic stimu- 
lation can enhance both short and long time correlation of heart beat. The ratio SD12 (points) 
in post stimulation state was higher than that in pre-state and the difference increased with lag 
number. The difference between the values of all three parameters before and after stimulation 
was found to be significant {p < 0.001). In order to find the relationship of these parameters 
with lag number M the method of Pade' approximant [23] was used. Assuming simple form of 
the Pade' approximant for SD's as 



a + bM 1 + 

y — = y (1) 

c + dM + ^ ^ 

the ratio of polynomial in M of degree one. Here Y = SDl, SD2 or SD12 and x — cl/c: /3 — b/a 

and ^ — d/c were taken as the new unknown parameters. The above equation was chosen 

by examining trend for small M and large M. For small M, Y increased hnearly and deviated 

at higher M. As shown below, the equ.(l) was found to be an excellent representation of the 

observed dependence of SD's on M (sohd hne in Fig.3). When expressed for small M the equ.(l) 

can be approximated as y = C + LM + QM"^ where L —x{P — l) and Q = - 7I/. It is to be 

noted that such variation of these parameters for small M were also found earher [18]. The 

values for L and Q were given in table 1. It was evident that the magnitude of slope and 

curvature of SDl after gyrosonic stimulation were increased considerably. On the other hand, 

the linear coefficient for SD2 increased slightly and curvature remained essentially unaltered. 

The ratio SD12 for different M for each subject was calculated and the mean value of the ratio 

was plotted with M in Fig.3 (points in lower curve). The data were excellently fitted by the 

equation (1) (solid curve) with the parameters value noted in table. 1. The ratio SD12 exhibited 

a larger change in both slope and curvature. The gyrosonic stimulation resulted respectively 

89% and 44 percent increase in the magnitude of curvature and slope of SD12. 



Table 1: The value of parameters X) 7 obtained from fit of eq.(l) with respective value of 
R^.The parameters L and Q are the coefiicient of linear and quadratic terms in expansion of Y 
in terms of M. Values of x, L and Q for SDl and SD2 are in second. 





a 


XxlO-2 


/3xlO-2 


7x10-2 


i?2 xlO-2 


L xlO-3 


Q xlO-^ 


SDl 


PRE 


1.34 ±0.02 


25.52 ±0.85 


2.15 ±0.11 


99.99 


3.14 ±0.13 


-0.67 ±0.02 




PST 


1.68 ±0.02 


27.25 ±0.96 


2.56 ±0.12 


99.98 


4.14 ±0.2 


-1.06 ±0.12 


SD2 


PRE 


3.54 ±0.03 


15.07 ±0.56 


3.12 ±0.17 


99.98 


4.23 ±0.31 


-1.32 ±0.1 




PST 


4.11 ±0.04 


13.74 ±0.53 


2.94 ±0.17 


99.97 


4.44 ± 0.34 


-1.31 ±0.11 


SD12 


PRE 


38.4 ± 0.24 


18.43 ±0.73 


8.74 ± 0.32 


99.97 


37.2 ±2.2 


-32.5 ±2.7 




PST 


40.04 ± 0.4 


24.58 ±1.1 


11.4 ±0.5 


99.95 


53.7 ±2.5 


-61.4 ±2.3 



Similar analysis was performed for individual subject and it was found that the equation 
(1) represented quite well with i?^ ~ 0.999. The values of L and Q for individual subject for 
pre and post-stimulation were depicted in Fig.4. Except for few (five to six) subjects the slope 
and magnitude of curvature increased after gyrosonic stimulation. Earlier study on subjects 
with cardiac illness had shown that the curvature of SDl and SDl 2 curves were much reduced 
compared to those for normal subjects [17]. Based on this result it can be argued that the 
augmentation of the curvature due to gyrosonic stimulation can then be taken as an indicator 



of better cardio-dynamics. The gyrosonic may be acting as a 'reward' signal for improvement 
of sympathovagal balance. 

The coefficient a and its distribution of detrended fluctuation analysis for subjects was plotted 
in Fig. 5a and 5b respectively for pre and post state of stimulation. The gyrosonic stimulation 
produced a decrease in the coefficient a for most of subjects and mean for the group is lowered. 
The mean of a for group changes from 0.93 ± 0.02 to 0.82 ± 0.02 after stimulation with a p 
value of 0.003. It is to be noted that no signiflcant changes were found in the above mentioned 
parameters in a group of subjects when we used various monotonous single frequency sound 
stimulation instead of gyrosonic stimuli. 

6 Conclusions 

The ECG data for short duration (10 min) is used to analyze the heart rate variability in time 
and frequency domain in order to assess the influence of novel gyro-sonic stimulation. It is 
found that the values of parameters SDl, SD2 and their ratio SD12 that quantify the Poincare' 
analysis of RRn interval become higher in post stimulation. The variation of these parameters 
with lag number is represented by an equation that fits excellently the data of SD's for group 
average and individual. The coefficients of linear and quadratic (curvature) term of SD12 
(and SDl) vs M relationship are enhanced due to gyrosonic stimulation. The equation (1) is 
an important findind as it provides the quantitative maesure of variation of SD12 with M - in 
particular the curvature of the plot which as appears from this and earlier [18] can be taken as 
a good measure of state of cardio-dynamics. The coefficient a of DFA is also decreased in post 
stimulation state. Gyrosonic stimulation of short duration reduces the heart rate and enhances 
SD and RMSSD. In frequency domain both LF and HF power are higher in post stimulation 
condition. The trend of the changes of most of the indicators of HRV supports the hypothesis 
that gyrosonic has the capacity to influence sympatho-vagal regulation in a favorable way 
and the gyrosonic can be considered as reward signal for improvement of autonomic activity. 
The importance of sympatho-vagal regulation on various diseases has already established. So, 
there is a possibility that the gyrosonic may act through improvement in autonomic regulation. 
However, more work on a larger group of patients with different health conditions are needed 
for establishing all aspects of gyrosonic. 
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Figure 1: Change of mean heart rate (HR),mean SD,mean RMSSD in RR-interval and mean 
Power of LF and HF before (PRE) and after (PST) gyrosonic stimulation. For all parameters 
the change is significant at p< 0.05. 




Figure 2: Poincare plot RRn+m vs RRjsi of one subject before (left) and after (right) gyrosonic 
stimulation for M = 1, 5, 9 and 13. 
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Figure 3: Variation of mean SDl and mean SD2 (upper) and mean SD12 (lower) with lag 
number M before (PRE) and after (PST) gyrosonic stimulation. 




Figure 4: Plot of slope L and curvature Q for individual subject before (PRE) and after (PST) 
gyrosonic stimulation. Unit of L and Q is in 10~^s. 




Figure 5: The DFA coefficient a of individual subject (left) and the probability distribution of 
a (right fig.) before (PRE) and after (PST) gyrosonic stimulation. 



